Transmission through a many-channel random waveguide with absorption 
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We compute the statistical distribution of the transmittance of a random waveguide with ab- 
sorption in the limit of many propagating channels. We consider the average and fluctuations of the 
conductance T = tvtH, where t is the transmission matrix, the density of transmission eigenvalues r 
(the eigenvalues of tH), and the distribution of the plane- wave transmittances T a and T a b- For weak 
absorption (length L smaller than the exponential absorption length £ a ), we compute moments of 
the distributions, while for strong absorption (L 3> £ a ), we can find the complete distributions. Our 
findings explain recent experiments on the transmittance of random waveguides by Stoytchev and 
Genack [Phys. Rev. Lett. 79, 309 (1997)]. 
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I. INTRODUCTION 



Interference between multiply scattered waves leads to 
strong fluctuations in the transmitted intensity through a 
disordered mesoscopic sampled While a theory of the pri- 
mary fluctuation phenomena was originally constructed 
in the context ftfj,electron transport through metals or 
semiconductors pel it was soon realized that the fluctu- 
ations are characteristic of systems involving multiple 
elastic scattering of any kind of waves, like sound, mi- 
crowaves, and light, and that they can described within 
the same theoretical framework.! Experiments with light 
or microwave radiation can be carried out with a very 
high accuracy, and allow for a precise verification of the 
theoretical predictions for the complete transmission jdis.- 
tribution, rather than of its mean and variance onlyETLl 

The relative importance of the fluctuations is deter- 
mined by the ratio N£/L, where N is the number of 
transverse propagating channels in the waveguide, L its 
length, and I is the elastic mean free path. To observe 
strong fluctuations, it is important to achieve as low val- 
ues of Nt/L as possible. This is difficult for optical ex- 
periments, because scatterers are weak ajid the typical 
number of transverse channels N is large B Moreover, in 
contrast to electronic systems, where the total flux is con- 
served, in optical systems the interference pattern may be 
affected as a result of loss or absorption. Absorption does 
not destroy the phase coherence of a transmitted wave; 
it merely rearranges the interference pattern, while the 
relative importance of the interference phenomena is un- 
affected. This is different from electronic systems, where 
inelastic scattering due to e.g. electron-electron interac- 
tions obscures the interference phenomena because of the 
presence of a large incoherent background signal. 

An important step in the pursuit of low values of 
N£/L for optical systems was recently reported by 
Stoytchev and Genack .0 They achieved N£/L sa 3 for 
microwave transmission through a copper tube with ran- 
domly placed polystyrene scatterers (a "random waveg- 







uide"). Apart from the length-dependence of the vari- 
ance of the transmittance, the measured transmittance 
distributioruatfrees surprisingly well with the theoretical 
predictions .eiO The agreement is surprising for a number 
of reasons. First, because the experiment is in a regime of 
strong absorption, the longest waveguide being approx- 
imately five times the exponential absorption length £ a , 
while the theory of Refs. p|-|Tl| was derived for waveguides 
without absorption. Second, because the variance of the 
transmittance depends sublinearly on L, which can not 
be explained within the existing theory for non-absorbing 
random waveguides. And third, because the experiment 
was compared to a theory for N£/L ^> 1, while Nl/L 
was not really large in the experiment. 

Several papers have dealt with the problem of trans- 
mission through absorbing— .waveguides with only one 
propagating mode (N — l).Lf~ll3 In order to analyze and 
explain optical or microwave experiments like that of Rcf. 

which are done in waveguides with many propagating 
channels, it is necessary that a theory of the transmis- 
sion fluctuations and the localization transition in many- 
channel random waveguides with absorption (N ^> 1) be 
developed. It is the aim of the present paper to present 
such a theory. 

The geometry of the random waveguide is quasi one- 
dimensional (width W much smaller than the length L). 
The relevant length scales are the elastic mean free path 
I, the exponential absorption length £ a , and the local- 
ization length £ — N£, where N ^S> 1 is the number of 
transverse channels in the waveguide, see Fig. |l|a. Al- 
though a priori the localization length £ is a property of 
the non-absorbing system, we find that it also governs 
the interference effects and the localization transition in 
a quasi one-dimensional random waveguide with absorp- 
tion. We assume that t <C £ a <C £, i.e. absorption is 
weak on the scale of a single scattering event, but it be- 
comes dominant before interference effects cause waves 
to localize. This is appropriate for experiments on the 
transmission of light and microwave radiation through 
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random waveguides. B~I3 The assumption £ a <C £ is cru ~ 
cial for our theory of the transmittance fluctuations in 
the localized regime L>^. 
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FIG. [j]. (a) The random waveguide that we consider here 
is quasi one-dimensional: its length L is much larger than its 
width W . The other relevant length scales are the wavelength 
A, the elastic mean free path £, the exponential decay length 
£ OJ and the localization length £ = N£, where N ^> 1 is the 
number of transverse channels in the waveguide. We assume 
A <C t -C Co ^ £• No assumption is made about the length L 
compared to £ a and £. (b) A thin slice of length 5L (left) is 
added to a random waveguide. 

Transmission through the waveguide is described by 
the N x N transmission matrix t, from which the three 
principal types of transmittances can be computed, 
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The transmittance T is the equivalent of the conductance 
for an electronic system. It is the transmitted intensity 
if the sample is illuminated through a diffusor (all chan- 
nels have equal incident power). The transmittances T a 
and T ao measure the total transmitted intensity and the 
intensity in channel b, respectively, if the sample is illu- 
minated through channel a only. (For N 3> 1 , this corre- 
sponds to plane- wave illumination.) The transmittance 
T a b is related to the speckle pattern, the configuration of 
randomly positioned dark and bright spots observed if a 
disordered sample is illuminated by a laser beam. The 
transmission matrix t can be decomposed into unitary 
matrices u and v and a matrix of transmission eigenval- 
ues r p (//=!,..., N), the eigenvalues of t't, 



t = u diag (tj 



1/2 



(2) 



Time-reversal symmetry implies that v — u T . For com- 
parison with the electronic case, we also address the case 
of broken time-reversal symmetry. The matrices u and v 



are uniformly distributed in the unitary group,0 as in the 
absence of absorptionJla To find the transmittance dis- 
tribution, it remains to find the statistical distribution 
of the transmission eigenvalues r M in the limit N 1 
corresponding to thick waveguides. 

This paper is organized as follows: In Sec. |l| we recall 
the scattering approach for the distribution of the tcaiis^. 
mittance, and its extension to absorbing systems ESEZI 
The transmittance distribution in the diffusive regime 
I <C L <C £ is considered in Sec. III. In IV we consider 



the crossover to the localized regime L > £. Sees. [II and 
[V primarily focus on the statistical distribution of the 
transmittance (or conductance) T. The distribution of 
the transmittances T a and T ao is discussed in Sec. [v| In 
Sec. VI, we discuss the relation ofr©ur work to the exper- 
iments of Stoytchev and Genack.0 We conclude in Sec. 



VII 



II. SCATTERING APPROACH 

The statistical distribution of the transmission ma- 
trix is obtained using a scattering approach similar the 
Fokker-Planck approach to the distribution of trans- 
mission eigenvalues in a disordered waveguide without 
absorption.li3~E2l The Fokker-Planck approach was also 
applied to the reflection -eigenvalues of a disordered ab- 
sorbing waveguide. ll3E3c3 Technical difhcultiesE-3 pre- 
vented a further generalization to the transmission eigen- 
vaiaes-|Of an absorbing waveguide beyond the case N = 
iJLata Starting from the random-matrix model of Ref. 
|l7| , we take a slightly different approach, which is ex- 
plained below. 

Consider a disordered waveguide with N propagating 
channels and length L, see Fig. |]. Its transmission and 
reflection properties arc described by the scattering ma- 
trix S. The 2N x 2N matrix S has the standard decom- 
position into N x N reflection and transmission matrices, 



S = 



r t' 
t r' 



(3) 



We now add a thin slice of width 5L to the waveguide 
and calculate the change of the transmission matrix t and 
the reflection matrix r. The slice has scattering matrix 
Si , which is parameterized aai3 



Si = 



n A 
h ri 



v' y/pv v' 'y/ru' 



(4) 



Here u, u', v, and v' are TV x N unitary matrices and 
t = diag (n, . . . , tjv) and p = diag {pi,...,p N ) are diag- 
onal matrices containing the reflection and transmission 
eigenvalues of the thin slice. In the presence of time- 
reversal symmetry ([3 = 1), one has v! = u T and v' = v T . 
A statistical ensemble of disordered waveguides is ob- 
tained by considering waveguides with different configu- 
rations of the scatterers. Following Ref. |l7|, we assume 
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that u, v! , v, and v' are uniformly distributed in the uni- 
tary group, and that the first moments of the diagonal 
matrices p and r are 

N~ 1 {trp) = SL/i, (5a) 
A^ 1 (tr r) = 1 - (SL/e + 5L/£ a ), (5b) 

where £ is the elastic mean free path and t a is the ballistic 
absorption length. The ballistic absorption length £ a is 
related to the exponential decay length £ a as £ a = 2£ 2 /£ 
Upon addition of the thin slice at the left end of the 
disordered waveguide, its transmission matrix t and re- 
flection matrix r are changed according to 

t -> i(l - r[r)~Hx, (6a) 
r n + t[r(l - r' x r)-Hi. (6b) 

The new transmission and reflection matrices t and r do 
not depend on t' and r'. Since we know the statistical 
distribution of the matrices t±, t'i, n, and r[ of the thin 
slice, and of the matrices t and r of the waveguide at 
length L, we thus can find the statistical distribution of 
the transmission and reflection matrices t and r of the 
waveguide at length L + 8L. In this way, one obtains a 
Fokker-Planck equation for the distribution of r and t. 

In non-absorbing random waveguides (£ a — > oo) 
this Fokker-Planck equation depends on the transmis- 
sion eigenvalues t\, . . . , tn only. It is known as.— tbs 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equationJlSE2l 
and it is one of-the major tools for the study of quan- 
tum transport S3 The DMPK equation has been gen- 
eralized to the reflection, eigenvalues of random waveg- 
uides with absorptionB-TEa For transmission through 
absorbing waveguides with many propagating channels 
(N ;§> 11, however, the Fokker-Planck approach proves 
useless :E£l The transmission eigenvalues do not decouple 
from the eigenvectors of tH and r'r, so that the num- 
ber of variables is of order A^ 2 , rather than TV. In this 
work, we take a different approach: we use Eqs. (g) and 
to derive a set of partial differential equations for the 
ensemble-averages of traces of (products of) r and t, with- 
out direct reference to the transmission eigenvalues r^. 
[We need to include the reflection matrix r, because the 
L-evolution of t depends on r, cf. Eq. (||).] A similar set 
of evolution equations for traces of the form tr (tH) n for 
the case of non-absorbing random waveguides has been 
derived from the DMPK equation in Refs. ^l] and from a 
microscopic theory in Ref. In the next two sections we 
present a detailed discussion of these evolution equations 
and their solution in the limit of random waveguides with 
many channels (N ^> 1). 

III. DIFFUSIVE REGIME 

In this section we consider evolution equations for 
traces of products of the reflection matrix r and the trans- 
mission matrix t, and obtain a solution as an expansion 



in \/N. We keep the ratio L/t fixed as we expand in 
1/N. Such a expansion is valid in the diffusive regime 
I <C L <C £, where £ = N£ is the localization length 
of the system in the absence of disorder. We explain the 
method by the computation of the averages (T) = (titH) 
and (ti r'r) to leading order in N, and then discuss the 
more general traces tr(vt) n , the density of transmission 
eigenvalues r, the weak-localization correction to the av- 
erage transmittance, and the transmittance fluctuations. 

A. Average of the transmittance T 

The simplest evolution equations are those for the av- 
erage transmittance (T) = (tr vt) and the average re- 
flectance (R) = (trrV). Combination of Eq. (||) and (||) 
yields 

d L {tvth) = - {r 1 +£a 1 ) (trth) +c £- 1 (trthtrr^r) 
+ S f3tl cp£~ 1 (trthr t r), (7a) 

d L (tr r^r) = -2(r 1 +^ 1 ) (trrV) + c^ -1 ((trr t r) 2 ) 
+ 7V + <5 /3il c^- 1 (tr(rtr) 2 ), (7b) 

where cp — P/(/3N + 2 — (3) and (3 = 1 (2) in the presence 
(absence) of time-reversal symmetry. Although Eq. (Q) 
does not form a closed set of equations from which the 
averages (tr tH) and (trrV) can be computed directly, 
it can be used to compute the transmittance distribu- 
tion in the limit of many channels (JV 3> 1). The rea- 
son is that the terms that couple the L-dependence of 
(tr tH) and (tr rV) to averages of traces with higher pow- 
ers of t and r are small by a factor of order N, so that 
their effect can be taken into account perturbatively. For 
non-absorbing waveguides, such large- A^ expansions have 
been studied in Ref. for the reflection properties of 
absorbing waveguides, a large-A^ solution according to 
these lines was given in Ref. ^8|. 

To find the leading large- N behavior of (tvtH) and 
(trrV), we retain only terms that are of order Af in the 
differential equation (0). A trace is counted as a factor 
N. To leading order in N, the average of, a product of 
traces equals the product of the averagesjH 2 ] while correc- 
tions are of relative order A^ -2 . Further, we may neglect 
(tr tHr^r) and (tr (rH) 2 ) with respect to (tr tH) (tr rH) 
and (trr^r) 2 . Then Eq. ([?]) simplifies to 

£d L (tr tH) = - (1 + 7- AT- X (tr r f r)) (tr tH) 

+ 0(1), (8a) 

£d L (trrH) =-2(1 +7) (trrH) + N + JV _1 (tr r f r) 2 
+ 0(1), (8b) 

where 7 = £/£ a - These equations are the same as those 
obtained from a diffusion-equation approach, neglecting 
the wavelike nature of the radiation. In the next subsec- 
tions, interference corrections will be taken into account 
by addition of the terms that we discarded as we sim- 
plified Eq. (@) to Eq. (§). The solution of Eq. (|) with 
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initial conditions (trtH) 
for L,£ a > i 

(T) = {trth): 



N, (trrV) 

i 

£ a sinh s ' 



at L = reads 

(9a) 
(9b) 



B. Traces of the form tr (tH) n and the density of 
transmission eigenvalues 



(R) = (tr r T r) = N - j- coths, 

where £ a = [ita/2] 1 / 2 and s = L/£ a . The length scale £ a 
is the classical exponential decay length for an absorb- 
ing random waveguide. In the weak absorption regime 
L -C £ a , Eq. (|) simplifies to Ohm's law (T) = N— (R) = 
£/L, while in the strong absorption regime L ^> £ a , the 
reflectance i? saturates at the value N — £/£ 0; while the 
transmittance T decays exponentially with decay length 



We now generalize the evolution equation (| 
trary traces of the form 



M x 



(tr xjxi . . 



■ x~ x r< 



to arbi- 



(11) 



(T) 



2£ 
— f 



(10) 



where the symbol can be t or r. These traces are im- 
portant for the density of transmission eigenvalues and 
for the distribution of the transmittances T a and T a b in 
the diffusive regime, see Sec. [v| Repeating the steps 
leading to Eq. M), we find 



£d,M x 



-(pi + l)(l + 7)Af sl 



k=l 



Xh M r ^ Xk+1 , 



fe=i 



,xj,_i-*'a;fc + i, 



E 



Pi + P?7 

iV 



k=2 

cyclic permutations 



0(1), 



(12) 



M x 



AO) 



3 = 1 



where pj = 1 (0) if Xj is r (t). The solution for L,£ a ^ £ and n = 2 reads 



<tr(^) 2 ) = f 
(tr^irV) = — 



^ /2s + coths 



4 sinh s 4 sinh s 
1 s coth s — 1 



4 sinh s 



4 sinh 3 



(tr(rV) 2 ) = N-j- 



3 coth s 



s 

coth s 



4 sinh s 4 sinh s 



(13a) 
(13b) 
(13c) 



r 



where as before s — L/£ a - For weak absorption, Eq. (IT 
agrees with results obtained from the DMPK equation 
while for strong absorption Eq. ( |l3| ) simplifies to 

(tr tVr) = (£/26,)e- L /*«, (14) 
(tr(rV) 2 ) = iV - (3£/2£ ). 

The averages (tr (tH) n ) correspond to moments of the 
density p(r) of transmission eigenvalues r. The density 
p(r) is recovered from the moments as the imaginary part 
of the Green function G(z) 

p(r) = n-hrG(T + iO), 

'tr (thy 



G(z) = (ix (z - tft)- 1 ) = 



E 

n=0 



(15a) 



For the calculation of p(r) we thus need the moments 
Af*,....t = N _1 tr (t't) n for all n. In principle this requires 
computation of all moments M xll ... lXm with m < n and 
with Xj being either r or t. As the number of possible 
moments M Xl ^_ >Xm proliferates exponentially fast with 
increasing m, this is not feasible. This is a fundamental 
difference with the set of moment equations for the case 
of a non-absorbing waveguide, where only moments of 
the form (tr (tH) m ) need to be taken into account. 



In the strong absorption regime L > („, however, a 
solution for p(r) can be found in closed form. Inspection 
of the general evolution equation ( |l2| ) for L ^> £ a shows 
that the leading behavior of tr (tH) n and tr (tH) n r^r has 
the asymptotic form 
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tr{{tk) n ) = a„^ n " 1 C" e "" L/eo + 
(tr (<ti)Vr) = Mi"" 1 ^-" 1 /*' + . 



(16a) 
(16b) 



where a„ and &„ are numerical coefficients. The dots 
indicate terms that are smaller by a large factor £ a /£ 
or exp(— L/^a). The evolution equation ( |l2| ) provides a 
recursion relation between the coefficients a n and b n , 



- J2 b ^ 



i 



(17) 



From Eq. (jlO|) we find a% = 2. The generating function 
F(z) — Xm=i a n z ™ 01 the coefficients a„ is the so-called 
"product-log" function (the principal value of the func- 
tional inverse of x — > xe x ), 



F(z) = ^a„z" = -2Plog (-z). 



(18) 



n=l 



The product-log function is re al on ly if its argument is 
larger than — 1/e. Using Eq. (15a), we find the Green 
function G(z), 



G{z) 



N 

z 



2£ 



Plog 



Le- L /t° 



(19) 



and hence the density of transmission eigenvalues p(r). 



My 
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FIG. ^. Density of transmission eigenvalues of a strongly 
absorbing random waveguide in the diffusive regime £ a >C 

The density of transmission eigenvalues is shown in 
Fig. [|. The existence of a maximum transmission eigen- 
value, r max = eLe~ L ^ a I '£ a <C 1 is quite different from 
the case of a non-absorbing random waveguide, where 
the support of pLz) extends from to 1 throughout the 
diffusive regime. E3 A common feature of absorbing and 
non-absorbing random waveguides is that the maximal 
transmission eigenvalue T max is a factor ~ L/l^> 1 larger 



than the average transmission eigenvalue. It is this broad 
support of the density p(r) that is responsible for many 
of the qualitative similarities of interference phenomena 
in strongly absorbing and non-absorbing systems. 



C. Weak localization correction 

Weak localization is a small negative-interference cor- 
rection to the classical transmittanceJ2j'E2l In electronic 
systems, the weak localization correction is suppressed by 
a time-reversal symmetry breaking magnetic field. In the 
framework of an expansion in 1/N, the weak localization 
correction is the 0(1) correction to an 0(N) average. 
Here we compute the weak-localization corrections ST 
and SR to the average transmittance (T) = (txtH) and 
reflectance (R) = (trr'r). 

Starting point is the exact evolution equation (Q), 
where we now keep all terms up to sub-leading order in 
N, 

£d L (trth) = - (l + 7 - A^trrV)) (tr i f i) 

- Sp,iN~ 2 {txr^r){txtU) 

+ 5 /3tl N- 1 (trthr^r) (20a) 
ld L (txr^r) = -2(1 + 7) (trrV) + 1 + N~ l {txr^r) 2 

- 8 0<1 N- 2 (txr^r) 2 

+ <WiV- 1 (tr(r t r) 2 ). (20b) 

We substitute the result of the previous subsection for 
(trt+tr+r) and (tr (r^r) 2 ) and solve Eq. @ for {txtH) 
and (trrV) up to order unity. As a result, we find 
the weak-localization corrections to the average reflec- 
tion and transmission, 



ST = 5 



0.1 



SR = S, 



0,i 



/ coth s — 2s 
\ 4 sinh s 
1 s coth s 

- -A = 

4 



4 sinh s 
- 1' 



(21a) 



(21b) 



4 sinh s 

The weak localization correction to the reflectance R 
has been calculated previously in Ref. ^8|. In the weak 
absorption regime, we recover the well known universal 
value SR = —ST = ^Sp^jEi while for strong absorption 
the weak localization correction reads 



ST 



-<3/3,i— e 



-L/C* 



SR 



1 



3/3,1- 



(22) 



D. Mesoscopic fluctuations 

Mesoscopic fluctuations of the transmittance T and the 
reflectance R are characterized by the variances varT, 
var R and the covariance cov (R, T) . Like the weak- 
localization cprrection, they are of order 1 in a large- 
N expansions Proceeding as in the case of the weak- 
localization correction, we find that varT, vari? and 
cov (R, T) obey 
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'9 L varT: 



ld L cov {R,T) 



idLV&r R 



-2(1 + 7- N^ 1 (tr rV» var T 

+ 2N- 1 (trth)cov (R, T) + 2iV~ 1 (tr (t t i) 2 r t r) 

+ 2^ (3il ^ 1 (trt t tr t t T t*r), 

-3(1 + 7- iV" _1 {tr rV)) cov (i?, T) 

+ iV" 1 (tr ^i) var i? - 4/3" 1 7V" 1 (tr ^trV) 

+ 4iV- 1 /3- 1 (tri t t(r t r) 2 ), 

-4(1 + 7- Af _1 (tr rV» var i? 

+ 4/3~ 1 iV- 1 (trr t r(l -rV) 2 ). 



(23a) 



(23b) 



(23c) 



The averages of the form (tr (tH) n (r^ r) m ) , can be computed from Eq. ([H 
we also need to know (titHrH T t*r), which satisfies (for 0=1) 



In the presence of time-reversal symmetry, 



£d L {tTthrh T t*r) = -4(1 + 7 - N~ l {tTr 1: r)){tTtHrH T t*r) 

+ N- 1 {trth)(2{tTth(r t r) 2 ) -4{trtHr^r) + (tri f t)) 
+ 3A^ 1 (trt t tr t r) 2 . 



(24) 



The solution of Eqs. 



2J) and (|2j) reads for L,£ a > 

2s 2 - 9s coth s + 12 
var T = - 



3s 2 



16 sinh 2 s 



16 sinh 4 , 



is coth s — 11 6s — 3s coth s — 3 



3s 2 



cov (i?, T) 



var R = — 



16 sinh s 
3 coth s 2s 2 coth s 



16 sinh s 
5s — 3 coth s 



8 sinh s 
3s 2 coth s 



3s 



16 sinhs 
1 1 
8 16 sinh 2 



16 sinh 3 s 
4s 2 - 3s coth s - 3 
16 sinh 4 s 

I 



16 sinh 5 s 



3s 2 



8 sinh s 



(25a) 
(25b) 
(25c) 



conductance fluctuations varT = —cov (R, T) = vari? = 
2/15/3, while for strong absorption Eq. ( |25| ) simplifies to 



0.2 




varT 
cov(i?, T) 
var i? 



L 2 
3 

~4/3 f 



-2L/£„ 



(26) 



FIG. |^. Fluctuations of the transmittance T in the diffusive 
regime L <C £. 

The L-dependence of the transmittance fluctuations, 
normalized to the average, is shown in Fig. ^. For weak 
absorption, L <C we find the well known universal 



The variance of R was obtained earlier in Ref. g8|. Note 
that the relative size of the transmittance fluctuations 
is equal in strongly absorbing and weakly absorbing sys- 
tems, varToc (L/£) 2 (T) 2 . 

The dependence of the fluctuations on the presence or 
absence of time-reversal symmetry is different from what 
one expects from the non-absorbing case and deserves 
some discussion. In non-absorbing disordered systems, 
the variance of the transmittance is decreased by a fac- 
tor 2 if time-reversal symmetry is broken, var T = 2/15/3. 
This universal 1/(3 dependence is well yuiderstood in 
terms of diagrammatic perturbation theorycra or random- 
matrix theory.E!] For a strongly absorbing system, in con- 
trast, the size of the transmittance fluctuations does not 
depend on the presence or absence of time-reversal sym- 
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metry, see Eq. ( p6| ) . This is remarkable, since the average 
transmittance in strongly absorbing waveguides exhibits 
the usual /3-dependent weak localization correction that 
is suppressed if time-reversal symmetry is broken. 

To see why the varT does not depend on the pres- 
ence or absence of time-reversal symmetry in absorb- 
ing systems, we consider the relevant pairs of the Feyn- 
man paths for the diffuson and cooperon contributions to 
varT. They are shown in Fig. |j. In a waveguide without 
absorption, both the diffuson and cooperon paths have 
equal weight. The cooperon is suppressed if time- reversal 
symmetry is broken, thus explaining the factor two reduc- 
tion of the fluctuations. In the presence of absorption, the 
diffuson and cooperon contributions no longer have equal 
weights, as paths for the cooperon contribution typically 
have a larger length. Therefore, unless the vertices 1 and 
2 are within a distance £ a , paths of the cooperon type are 
suppressed. Since the diffuson paths of Fig. |]a can have 
vertices arbitrarily far apart, the cooperon contribution 
to varT is a factor £ a /L smaller than the diffuson one. 
This explains why the breaking of time-reversal symme- 
try has no effect on the transmittance fluctuations. The 
situation is different for the weak-localization correction 
to the average transmittance. The appropriate Feynman 
paths, which are shown in Fig. ||, contain a closed loop 
with a typical length of order £ <§; l a . As a result the 
weak-localization correction to the average transmittance 
(T) is not suppressed by absorption. 




(a) (b) 

FIG. Feynman paths for the diffuson (a) and cooperon 
(b) contributions of the variance of the transmittance T . 
There is only a cooperon contribution to var T if the vertices 
1 and 2 are within a distance £ a of each other. For diffusons, 
the vertices can be arbitrarily far apart. Hence for i > 
only the diffuson contribution to var T survives. 




FIG. |5| Two interfering Feynman paths for the weak- 
localization correction to the average transmittance T. As 
the typical length of the closed loop is of order f < 4, the 
weak-localization correction is not suppressed by absorption. 



E. Breakdown of the large-TV expansion 

Throughout this section we have assumed that interfer- 
ence corrections are small, so that we can treat them per- 
turbatively. This assumption is bound to break down as 
the length L of the waveguide increases, since the relative 
size of the interference corrections to transmission prop- 
erties increases with L. For the transmittance T, both 
the fluctuations and the weak-localization correction be- 
come comparable to the average when the length of the 
waveguide approaches the localization length £ = N£, cf. 
Eqs. ( p2[ ) and (p6[). One verifies that higher order cor- 
rections in a 1/N expansion become comparable to (T) 
as well as L — > £. This breakdown of perturbation the- 
ory has the same origin as the corresponding breakdown 
of perturbation theory in systems without absorption: it 
signals the onset of localization. A theory of the trans- 
mittance distribution in the regime L > £ where the per- 
turbation theory of this section is not valid, is presented 
in the next section. 



IV. CROSSOVER TO LOCALIZED REGIME 



The large- N perturbation theory of the previous sec- 
tion breaks down as the length L approaches the local- 
ization length £. For non-absorbing random waveguides, 
a theory of the crossover from the diffusive regime into 
the localized regime requires a true technical tour dc 
force,Ej~0 because of the intrinsically non-perturbative 
nature of the crossover. Although the situation for ran- 
dom waveguides with absorption looks similar — it has 
the same divergence of perturbation theory as L ~ £ 
— , it is not. The reason is the existence of the small 
parameter or l/"fN 2 . (Note that the ratio £ a /£ 

depends on the absorption properties of the waveguide 
and its width, but not on its length.) In this section, we 
use the smallness of £ a /£ to compute the distribution of 
the transmittance for lengths L comparable to, or greater 
than the localization length £. 

A crucial distinction between the case of absorbing 
and non-absorbing waveguides is the saturation of the 
reflectance distribution for absorbing waveguides for L > 
£ Q . Although the relative size of the transmittance fluctu- 
ations keeps growing as L exceeds the absorption length 
£ Q — eventually causing the breakdown of perturbation 
theory — , the correlations between reflection and trans- 
mission properties saturate. In particular, we find that 

• Correlations between reflection and transmission 
are smaller than the product of the averages by 
a factor ~ S,a/N£_, irrespective of the length L. 
Hence they can be neglected if £ a <C £. (For non- 
absorbing waveguides, correlations are smaller by a 
factor L/N£, and can not be ignored for L ~ £.) 
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l r are a factor 
or tr (t+t)"tr (thy 



• Traces tr (tH) n r\t T t*) 
smaller than tr (th) n+m 

A derivation of these properties from the appropriate evo- 
lution equations is straightforward. We find it more in- 
structive, however, to give an argument in terms of Feyn 
man paths. 



2- - ■ 





FIG. |6|. Feynman paths contributing to correlations be- 
tween the reflectance R (solid path) and the transmittance 
T (dashed path). The solid path, which contributes to the 
reflectance R, does not penetrate the sample further than the 
decay length £ a . Therefore, correlations between R and T 
saturate at L > £ a . 

Feynman paths giving rise to correlations between R 
and T are shown in Fig. ^. The solid path, which con- 
tributes to R, does not penetrate the sample more than 
a decay length £ a . As a consequence, the relative size 
of the correlator cov(R,T) saturates at L ~ £ a . At this 
length scale, perturbation theory is still valid, and we 
find cov(i?,T) - (£ a /N£){R){T) < (R)(T) irrespective 
of length. The same argument applies to arbitrary cor- 
relators involving tr (r'r) n , since their relative sizes sat- 
urate at L ~ £ Q as well, and to the traces of the form 
tr {tH) n r^(t T t*) m r. 

The formal framework in which the crossover to the 
localized regime is described, is the so-called "thick-wire 
limit" , in which the limit N — > oo is taken, keeping 
<j = L/£ fixed. In the previous section, the ratio L/i 
was kept fixed, rather than L/£. The thick-wire limit is 
used in the field-theoreticj-description of localization of 
in disordered waveguides. E3t3 Note that the thick-wire 
limit is unable to address length scales L < £ a corre- 
sponding to the onset of absorption, since it corresponds 
to a = — * if N — > oo. In order to remove the 
classical exponential decay of the transmittance due to 
absorption, we consider traces of the form 



Fn 



(27) 



The differential equations for the cr-dependence of the 
averages of the traces F n are derived from the exact evo- 
lution equations for the averages of (products of) traces 
tr (tH) n (r'r) m , which are constructed as discussed in Sec. 
EJ. As an example, we consider the evolution equation 
for (Fi) in the thick- wire limit. Hereto we rewrite the 
evolution equation (^) for (tr tH) in terms of a = L/£ 
and Fi = (2£/£ a ) exp(— L/£ a )tr tH and take the limit 
N — > oo. The resulting equation reads 



5 CT (F 1 ) = ^, 1 ((F 1 )-|(F 1 )) ) 



(28) 



where F[ — (2£/£ a )e L '£ a tTtHr^r. The evolution equa- 
tion for (F[) in the thick-wire limit takes a particularly 
simple form, 



(U0d a {Fl) = -2{F{) + (l/2){F 1 



(29) 



Here we used that (tr (1-rV) 2 ) = £/2£ Q for L > £ a . The 
left-hand side of Eq. (|29|) vanishes, so that (F[) = (Fi)/4. 
Substitution into Eq. (|28|) yields 



da(F 1 



(30) 



We use Eq. (|f0| ) for the initial condition, (F\) = 1 for 
er — > 0. (Notice that this is not the ballistic initial con- 
dition. Because of the order of limits taken, the initial 
condition a — > has to be evaluated within the strong 
absorption regime.) Hence 



(F 1 ) = e X p{-L6p, 1 /20. 



(31) 



For ( a « L « ( this agrees with the weak localization 
correction to (T) in the diffusive regime for strong ab- 
sorption, cf. Eq. (|H). 

Let us now consider the L-dependence of a general av- 
erage (lXj=i Fnj) m the thick-wire limit. Following the 



derivation of Eq. (30), we find the general evolution equa- 
tion 



MllM =EEi f <ViII f ' 



\3=i 



W =1 / a^O 



3=1 i=l 



fc=i 



i<3 \ 



fe=l 



3 = 1 \k=l 



3 = 1 



(32a) 
(32b) 
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The set of linear differential equations (32a) is readily 
solved. Its solution reads for JT n>j < 3 



/ J? \ 


lt r/2 


(ooa) 


(F?) 


= cosh (a/2)e- 5 ^ 1 ' 7 1 


(33b) 


(F2) 


= sinh (cr/2)e-^' 1 ' T , 


(33c) 


(F?) 


= | [cosh (3(t/2) + 2] e- 35 ^ 2 , 


(33d) 


(FiF 2 ) 


= isinh(3(j/2)e- 3,5 ' j ' icr / 2 


(33e) 


(Fs) 


= i [cosh(3cr/2) - l}e- 35 ^ 1,j/2 . 


(33f) 



the waveguide is illuminated through a diffusor. Using 
the results of the previous two sections, it is only a small 
step to find the distributions of the transmittances T a 
and Tab, which describe a situation in which the waveg- 
uide is illuminated through a single channel only (e.g. 
by a plane wave). In fact, it is sufficient to consider the 
transmittance T a , since._£he distribution of T a t follows 
directly from that of T a 



P{T ab ) = N 



dT a , s / NT ab 
yip(T a )expf-^ 



(38) 



Deep in the localized regime it>1, the solution of Eq. 
(|32a|) has the asymptotic form 



log(F") = \a [n(n - 1) - 2nS p . 1 ] - log. 



(34) 



In the diffusive regime L -C £, the fluctuations of the 
transmittance T are much smaller than the average. The 
variance of the transmittance T in the strong absorption 
regime can easily obtained from Eq. ( |33|) and agrees with 
the result of Sec. [II. We also note that the third cumu- 



lant (T 3 ) c behaves as 

{T 3 ) c /{T) 3 = UL/i) A + 0{L/£,f 



(35a) 



which is qualitatively different from the case of a non- 
absorbing waveguide, where ZX 3 !,/^) 3 ~ (L/£) 5 for 
f3 = 1 and - (L/£) 6 for (3 = 2£§E3 

For L ;§> £, the fluctuations of T are much larger than 
the average. Therefore, the average and variance of the 
transmittance are no longer sufficient to characterize the 
distribution. Instead, the transmission distribution is 
log-normal, like for a waveguide without absorption. This 
can be understood from the fact that for a log-normal 
distribution of T, the nth moment (T n ) reads 



log(T") = n(logT) + (n 2 /2)var \ og T. 



(36) 



Using the asymptotic result (|34j) , together with Eq. < \27\j 
we see that Eq. (5q) is satisfied for all n with 



(logT) 
var logT 



-£/£« 

L/2£. 



(L/40(l + 2^ ;1 ), 



(37a) 
(37b) 



[We neglected the— second L-independent term on the 
r.h.s. of Eq. (B^).E3] For a waveguide without absorp- 
tion, on the contrary, the average and variance of log T 
are approximately equal and depend on the localization 
length £ only: (logT) = -2L/(3£ and var logT = 4L//3£. 
The insensitivity to (3 of the transmittance fluctuations 
for a waveguide with absorption was already discussed in 
Sec. fn. 



V. TRANSMITTANCES T A AND Tab 

So far we have mainly considered the statistical dis- 
tribution of the transmittance T, which is appropriate if 



The statistical distribution of the transmittances T a 
and T a b for non-absorbing random waveguides has been 
calculated in Refs. 0-|ll|. For L <C £, the distribution 
of T a is Gaussian with mean (T a ) = N~ X (T) = l/L, 
variance varT Q = 2£/3NL, and non-Gaussian tails. The 
non-Gaussian features of the distribution become more 
pronounced as the length of the waveguide approaches 
the localization length £ = N£. In the localized regime 
has a log-normal distribution, with the same 
mean and variance as the transmittance Tr3 

Let us now discuss the effect of absorption on the dis- 
tribution of T a . (The variance of T a in the presence of 
absorption has also been considered in Refs. [h] and 
In terms of the unitary matrix u that diagonalizes ttf [cf. 
Eq. (|2j)] and the transmission eigenvalues r M , the trans- 
mittance T a reads 



(39) 



The matrix u is uniformly distributed in the unitary 
group. For large N, we may consider the matrix elements 
u aix as independently distributed Gaussian random num- 
bers with zero mean and variance (|u a /j| 2 ) = l/N. Iixlhis. 
limit, the Laplace transform F(z) of P{T a ) becomealJ'EiJ 



F{z) 



dT n e 



-NzT a 



P{T a ) 



(40) 



In the strong absorption regime L » („, we can sub- 
stitute the results of Sec. IV for the averages of moments 
of the transmission eigenvalues. The result is 



F(z) = J2 zn exp 



n=0 



L L L 

n(n-l)--n6p, 1 --n- 



(41) 



Transforming back, we find that T a is log-normally dis- 
tributed with mean and variance 

(logT Q > = - (1 + 2Sp,i)^ - logiV, (42a) 



var logT a = — . 



(42b) 
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The origin of the remarkable simple log-normal distri- 
bution of T a in the strong abspeption regime is under- 
stood from a simple argument x3 In the strong absorp- 
tion regime, diffusion takes place only inside a volume 
£ Q -C L, and motion is quasi-ballistic on larger length 
scales. Inside a segment of £ a , absorption plays no role, 
and the distribution of the transmittance T a is Gaussian, 
with mean ~ £/£, a and variance ~ ^ 2 /££a, as in the case 
of a non-absorbing random waveguide. Since the dynam- 
ics on larger length scales is quasi-ballistic, the transmit- 
tance T a of a waveguide of length L is the product of 
the individual transmittances of uncorrelated segments 
of size £ a . From the central limit theorem we then imme- 
diately derive that for L ^> £ a , the transmittance T a is 
log-normally distributed with variance var logT a ~ L/£. 



where s = L/£ a . Eq. (E^) is valid for both weak 
and strong absorption. In the weak absorption regime 
L <C £a , Eq. reduces to the known results for a non- 
absorbing waveguideoEJ 

9 T 1 R T 2 

varT a = ^(T a ) 2 , (T 3 ) c = i^<T a ) 3 . (44) 

whereas for strong absorption we find 

T "\T, 2 
varT a = -(T Q } 2 , <T 3 ) C = — <T a ) 3 . (45) 

The strong absorption limit of var T a was previously ob- 
tained in Refs. 




12 3 4 
T /(T ) 

FIG. [?]. Distribution of the transmittance T a for L/£ = 0.1, 
0.5, 2, and 10 in a strongly absorbing random waveguide. In- 
set: The same, but on a logarithmic scale. 

The distribution of T a in the strong absorption regime 
is shown in Fig. ^ for several values of £/£. In the local- 
ized regime L 3> £, the distribution of iVT a is the same as 
that of the transmittance (or conductance) T. However, 
the log-normal distribution ( fi"^ ) of the transmittance T a 
is valid for all L 3> £ a , he. before the onset of localization. 

In the diffusive regime ^ < L < the distribution of 
T a is sharply peaked around the average (T a ) = A r_1 (T). 
Fluctuations around the average are described by the cu- 
mulants 

varT a = ^-AT a ) 2 [2s + coths - -^ T -) , (43a) 
44 \ sinh s J 

6s 2 + 5scoths — 1\ , . 

r-2 , (43b) 

sinn s / 
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FIG. [|. Distribution of the transmittance Tab f° r L/£ = 
0.1, 0.5, 2, and 10, in a strongly absorbing random waveguide. 
The limiting negative exponential distribution for L -C £ is 
shown dotted. 

The distribution of the transmittance T a b in the strong 
absorption regime is easily obtained using Eq. (^3q). It 
is shown in Fig. |^ for several values of L/£. For L <C £, 
the distribution of T a b is negative exponential (Rayleigh) . 
In the localized regime L 3> £, the distribution is log- 
normal. 



VI. EXPERIMENTS 

In this section, we discuss the application of our work 
to the recent microwave experiments by Stoytchev and 
Genack.Lj In these experiments, the distribution of the 
transmittance T a was measured in a cylindrical copper 
tube with randomly placed polystyrene scatterers. The 
presence of many scatterers in a narrow waveguide results 
in a relatively small localization length (£ ~ 5m, while 
the length L of the waveguide varies between 50cm and 
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2m), but at the cost of strong absorption (£ a w 30cm). 
Lacking a theory for absorbing random waveguides, the 
authors of Ref. [?] considered the distribution of the trans- 
mittance divided by its mean, T a — T a /(T a ), and com- 
pared it to jthe_.theoretical predictions for non-absorbing 
waveguides n~i3 Because of the absorption, the average 
transmittance (T a ) can not be used to determine the ra- 
tio L/£. Instead, in Ref. 0, L/£ was computed from the 
formula var 7^ = (2/3) (L/£), which is valid for weakly 
absorbing waveguides only [cf. Eq. (Q)]. 

Surprisingly, the measured probability distribution of 
T a was found to follow the theoretical predictions for non- 
absorbing waveguides in the diffusive regime quite accu- 
rately, provided one uses the ratio L/£ obtained from 
var 7^, as explained above. This is surprising, because 
the waveguides of Ref. || are strongly absorbing, L/£ a 
ranging from 2 to 6, and because the experiment is close 
to the localization threshold (£,/L up to 3). The agree- 
ment is not so good if L/£ is directly extracted from the 
experimental parameters. In particular, the measured 
variance of T a was found to depend sublinearly on L, 
while the theory for a non-absorbing waveguide predicts 
a linear .L-dependence (or a superlinear L-dependence if 
localization effects are taken into account). 
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FIG. § (a) The ratio C(L) = <T a 3 ) c /(varT a ) 2 versus L/£ a 
in the crossover between weak and strong absorption, (b) The 
ratio (£/L)vaxT a versus L/£ a . The errorbars (a) and squares 
(b) indicate the results of the experiment of Ref. 0. 

How do these observations compare to our results 
for the transmittance distribution of random waveguides 
with absorption? Following Ref. [?], we use the ratio 
C(L) — (7^ 3 ) c /(var7^) 2 to characterize the distribution. 
The theoretical predictions for C(L) and (£/L)var7^ as 
a function of L/£ a in the diffusive regime are shown in 
Fig. | [see also Eq. ©]. We find that 

• as a result of absorption, var 7^ crosses over from 



§L/£ to \LH for 1 < L/^a < 10, 

• the ratio C(L) crosses over from C(L) = 12/5 for 
weak absorption to 3 for strong absorption, the 
crossover occurring for 3 < L/£ a < 30. 

The length L of the random waveguides in the experi- 
ment lies between 2£ a and 6£ a . Hence, the waveguides 
are long enough to fully observe the sublinear behavior 
of var T a , but too short to see a significant enhancement 
of C(L) above the weak absorption limit C(L) = 12/5. 
(Note that the experiment indeed shows a slight enhance- 
ment of C(L) for the longest waveguide.LJ) We expect 
that higher cumulants oiT a , if properly normalized, show 
the same "postponed" crossover behavior as {T^) c - This 
could explain why the entire transmittance distribution 
of Ref. [?] agrees with weak-absorption theory rather than 
with strong-absorption theory. 

While the above considerations offer a qualitative ex- 
planation for the most striking experimental observa- 
tions, the theory presented here fails to account for the 
measurement of var T a versus L quantitatively. This is 
clearly illustrated by comparison of the experimental re- 
sults of Ref. fj] and the theory in Fig. |^b. Notice that 
the fact that the theoretical curves in the figure were de- 
rived for L <C £, while the ratio £/L is not necessarily 
large in the experiment (£/L down to 3), cannot explain 
the difference, as localization effects increase fluctuations 
rather than decrease them. For the longest waveguide 
(L/£ a w 6), Stoytchev and Genack find var T a w 0.43L/£, 
whereas we find var T a > \L/^ for all L. The reason for 
this discrepancy is not known. A quantitative compari- 
son with the experiment for C(L) is difficult because of 
the uncertainty of the experimental data (see Fig. ||a). 

According to Eq. ( ^3| ) or Fig. |^a, absorption causes 
the distribution of the transmittance to deviate signif- 
icantly from that of a non-absorbing waveguide only if 
L/£ a <^ 10. It should not be difficult to verify this ex- 
perimentally, e.g. by a decrease of £ a due to the addition 
of strongly absorbing scatterers to the waveguide. We 
also find that the size of the fluctuations depends lin- 
early on L/£, both in the weak and strong absorption 
regimes, but with different slopes (2/3 and 1/2, respec- 
tively). While the experiment of Ref. [7] confirms the ini- 
tial -L-dependence with slope 2/3 as well as the deviation 
for L ~ £ a , more experimental input is required to ver- 
ify the linear dependence with slope 1/2 in the strongly 
absorbing regime. 

VII. CONCLUSION 

We have computed the statistical properties of the 
transmittance of a multichannel random waveguide with 
absorption. We have considered both the diffusive 
regime, where the transmittance distribution is Gaussian 
and the localized regime, where the transmittance distri- 
bution is log-normal. Our main findings are summarized 
in Table | 
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no / weak 
absorption 


strong 
absorption 


diffusive 


1 < L < £_ a 


i a < l < e 


P(T) 
(T) 
<5T 
varT 
P(r) 


Gaussian 
bimodal 


Gaussian 

(2C/e a ) e - L/ «» 

W =eL{T)/2i 


localized 






P(T) 
(logT) 
var log T 


log-normal 

-2L//3£ 
4L//?£ 


log-normal 
-£/£> + ©(£/£) 
£/2£ 



ACKNOWLEDGMENTS 



We thank C. W. J. Beenakker, T. Sh. Misirpashaev, 
and N. M. Shnerb for stimulating discussions. We thank 
A. Z. Genack and M. Stoytchev for correspondence and 
communication of the experimental data. This work was 
supported by the "Stichting voor Fundamenteel Onder- 
zoek der Materie" (FOM) and by the "Nederlandse or- 
ganisatie voor Wetenschappelijk Onderzoek" (NWO) and 
by the NSF under grants no. DMR 94-16910, DMR 96- 
30064, and DMR 94-17047. 



TABLE [ij Statistical properties of the transmittance T for 
weak and strong absorption. The relevant length scales are 
the the localization length £ = N£ and the exponential decay 
length £ a = [lla/2] 1 / 2 , where N is the number of channels, 
£ the elastic mean free path, and £ a the ballistic absorption 
length. 



To present, no optical or microwave experiments have 
been able to address the localized regime L > £, where 
the fluctuations of the transmittance are larger than the 
average. The mere presence of strong absorption modi- 
fies the localization transition, but does not suppress the 
fluctuations. As stronger scatterers tend to be stronger 
absorbers, it might be easier from the experimental point 
of view to construct a random waveguide with L > £ in 
the strong absorption regime than a comparable random 
waveguide without absorption. 

Our work assumes that absorption sets in before local- 
ization, i.e. that £ Q -C £. The assumption was necessary, 
because we use £ Q /£ as a small parameter. It is appro- 
priate for the experiment of Stoytchev and Genack, in 
which the ratio £ a /£ is less than 0.1. A theory of the lo- 
calization transition in random waveguides that is non- 
perturbative in £ a /£.— which would require a non-linear 
cr-model formulation^ remains to be developed. 

To conclude, we would like to remark that the prob- 
lem of transmission through random media in the pres- 
ence of absorption is akorinrportant for the problem of 
"directed localization" .e3e3 Directed localization refers 
to the localization transition in a disordered ring with an 
imaginary vector potential, first introduced by Hatano 
and Nelson.c3 This model is relevant for e.g. the pinning 
of vortices to columnar defects in superconductors, or 
for problems in population biologyO In a purely one- 
dimensional geometry (N = 1), the support of the spec- 
trum in the complex plane is entirely determined by the 
transmittance of the disordered, ring at complex energies, 
i.e. with absorption (or gain).cll Such a relation may also 
exist for the multichannel waveguides considered in this 
paper. 
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